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Abstract. A noisy Gaussian channel is defined as a channel in which an input field mode is 
subjected to random Gaussian displacements in phase space. We introduce the quantum fidelity 
of a Gaussian channel for pure and mixed input states, and we derive a universal scaling law of 
the fidelity for pure initial states. We also find the maximum fidelity of a Gaussian channel over 
all input states. Quantum cloning and continuous-variable teleportation are presented as physical 
examples of Gaussian channels to which the fidelity results can be applied. 



1. Introduction 

A principal aim of quantum information theory ^ is to determine the ultimate 
limits on capacity, entropy, or fidelity of quantum information transmitted in noisy 
quantum channels [21 EH HI El El • Quantum channels with noise are examples 
of open quantum systems that interact with an environment. The environment 
produces classical and/or quantum fluctuations that generally degrade the input 
states. The channel is described by a linear map p t—. > Q(p), which takes an input 
state p to an output state &(p)- The efficacy of the channel can be characterized by 
a quantum fidelity J-(p, 7), which measures the quality of transmitted information 
as a function of the input state p and the noise parameter 7 that describes the 
environment. 

A special class of noisy quantum channels consists of bosonic channels with 
excess noise described by random Gaussian shifts in phase space [HI [HI El ITT] . 
Such quantum Gaussian channels have attracted considerable attention in the the 
framework of quantum information with continuous variables. Recent investiga- 
tions devoted to bosonic Gaussian channels have addressed the following problems: 
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lower bounds on the capacities ^U], the question of multiplicativity of the maximal 
p-norm output purities ^21) an d the Holevo- Werner additivity of entropies at the 
output of the channels [T3] . 

The purpose of this paper is to describe the quantum fidelity of Gaussian chan- 
nels, using as a tool phase-space methods related to Wigner and Weyl functions. A 
general scaling law for quantum fidelity for pure input states is derived. We argue 
that this scaling law reflects the duality of the Wigner and Weyl phase-space dis- 
tributions related to sub-Planck 14 and large scales of the quantum state p (for 
an extended discussion of this duality see Ref . ^3] ) . We show that if the devices 
processing quantum information can be built from linear optical elements in an 
extended Hilbert space, the resulting Gaussian channels have important physical 
applications. We discuss quantum cloning and a continuous- variable teleportation 
protocol as examples of Gaussian noise channels. 

The paper is structured as follows. In Sec. ED the description of quantum 
channels in terms completely positive maps (CPM) is recalled. In Sec. ED the 
Gaussian channel is defined in terms of random Gaussian shifts and is related to 
Wigner and Weyl functions in phase space. In Sec. E] the quantum fidelity of 
a Gaussian channel is introduced for pure and mixed input states. In Sec. ED a 
universal scaling law of the fidelity for pure initial states is derived, and in SecFOwe 
find the maximum fidelity of a Gaussian channel over all input states. Section 171 
presents as examples explicit calculations of the fidelity for number states and 
squeezed states. Section |HD applies two different kinds of fidelities for mixed input 
states to a thermal input. In Sec. ED we discuss quantum cloning and continuous- 
variable teleportation as physical examples of bosonic Gaussian channels. Some 
final thoughts are presented in Sec. 110. 1 



Consider a quantum system described by a Hilbert space TC, with a given 
density operator p. A noisy quantum channel is a linear and trace-preserving map 
of the quantum state p, given by the Kraus decomposition 



where the completeness condition ]TV KjKi = I makes the map trace-preserving 
[T5| . We call p the input state and the transformed &(p) the output state. 

The overall system-environment state is described by a density operator that 
evolves unitarily through the quantum channel. Most often, one is interested in 
the system alone, which is described by a reduced density operator obtained by 
tracing over the environment degrees of freedom. As a result of such reduction, 
a quantum channel with noise is characterized by a linear, trace-preserving, and 



2. Quantum channels and completely positive maps 
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completely positive map. A completely positive map is defined in the following 
way: if the system undergoes the dynamics described by and a reference system 
R of arbitrary Hilbert-space dimension experiences no dynamics, a situation de- 
scribed by the overall superoperator $(8>Xr, where Ir is the identity superoperator 
for the reference system, then an arbitrary (potentially entangled) joint state of 
the system and reference system is mapped to a positive output state; i.e., the 
superoperator $ <8> Tr maps positive operators to positive operators. 

A completely positive map is required to describe reduced dynamics because 
it implies that such dynamics arises from a unitary evolution of the system and 
an environment, 

= Tr^[l7(p <8» |e><e|)C^] . (2) 

Here the environment degrees of freedom are denoted by E and \e) is some initial 
state of the environment. The following statements about a linear and trace- 
preserving map are equivalent: (i) <I> is completely positive, (ii) $ has a Kraus 
decomposition as in Eq. (fT|). and (iii) $ is the reduced dynamics for some system- 
environment unitary dynamics as in Eq. (J2J). 



3. Bosonic Gaussian Channels 

What we mean by a bosonic Gaussian channel in this paper is a completely 
positive map that results from zero-mean random Gaussian shifts in the phase 
space of a bosonic mode [HJ. This CPM acts on input states p in the following 
way: 

<&{p) = j d 2 aQ{a)D(a) P D\a). (3) 

Here D(a) = exp(aa^ — a*a) is the displacement operator in the phase space of 
the bosonic mode, where the bosonic creation and annihilation operators obey the 
commutation relation [a, a'] = 1, and the classical-noise Gaussian distribution, 

g(a) = —e~ 2 ^ 2 /\ (4) 

has zero mean and phase-space variance a 2 = 7/2. The complex number a = 
(q + ip) /y/2 corresponds to the phase-space points of a single-mode harmonic 
oscillator described by position q and momentum p. Notice that the Gaussian 
channel of Eq. (j3J) is a unital map; i.e., it leaves the identity operator unaffected. 
More importantly, this Gaussian channel preserves the expectation value of a; i.e., 
this kind of Gaussian channel cannot have any overall attenuation or amplication. 
The Gaussian channel © is in the Kraus form if we make the identification 
Ki -» ^/g(a)D(a). 

In Sec. E0 we show that such Gaussian channels can be implemented in the 
framework of linear optical elements in an extended Hilbert space. In such cases 
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the CPM © can be derived from the general relation We will discuss quan- 
tum cloning and continous-variable teleportation as physical examples of bosonic 
Gaussian channels resulting from a reduction of quantum systems that interact 
with an environment. 

In our discussion of the Gaussian channel, we will analyze various phase-space 
overlaps. We use the following definitions of the phase-space Wigner and Weyl 
functions. The Wigner function of an arbitrary density operator p is IT given by 

W p (a) = -Tr[ pD(a)(-lf a D\a)] = [ e aP *~ a ^ C p (f3) , (5) 

7T L J J TT Z 

where C p (a) = Tr[pD(a)] is the Weyl characteristic function and p is an arbitrary 
density operator. We find it useful to rewrite the Gaussian CPM Q in terms of 
the phase-space Wigner functions of the input and the output states: 



W Hp) {a) = J d 2 PG(a - (3) W p {(3) . (6) 

The counterpart of this relation in terms of the Weyl functions of the input and 
output states is 

C Hp) (a)=e-^ 2 / 2 C p . (7) 



4. Fidelity of a Gaussian channel 

4.1. Fidelity for pure input states 

If the input state to the Gaussian channel is a pure state p = \^)(^\, the 
channel fidelity is given as a quantum overlap between the input and the output 
state. In this case the channel fidelity is 

^(tt l7 ) = <*|*(p)|*>. (8) 



Using the phase-space Wigner functions of the input and output state, we can 
rewrite this fidelity phase-space overlap 

(a) = it d 2 ad 2 PG{a- (3)Wy{a)W*(P). (9) 

Another useful form for this fidelity, this time in terms of the Weyl functions, 
comes from substituting the output state © directly into Eq. |H|): 



^(¥,7) = / d 2 ag{a)\Cm{a)\ 2 



(10) 
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4.2. Fidelity for mixed input states 

An appropriate measure for assessing the fidelity of an mixed input state is the 
entanglement fidelity JE]; which is defined in the following way. Imagine that the 
input mixed state p is purified to a state of the original mode and a reference 
mode R. An example of such a purification is 

W) = Vp®IrY1 i n ) ®l n ) = ^2VP~\ n ) ®\ n ) ■ ( n ) 

n n 

The input mode is now entangled with the reference mode, and the purified state 
belongs to an enlarged Hilbert space TL ® Hr- We now suppose that the origi- 
nal mode is subjected to the Gaussian channel while the reference mode is left 
untouched. The resulting output state is 

$®Jji(MM) = J d 2 aQ{a)D{a)®I R \iP){iP\I R ®D\a). (12) 

The entanglement fidelity is now defined to be the fidelity of this joint output state 
with the purified input state: 

^(p,7) = <V|*®Xfl(|^)(^|)^>. (13) 

The entanglement fidelity is a property of the system state p even though the 
purification is not unique. The entanglement fidelity reduces to the fidelity of 
Eq. (jHJ) in the case of pure input states. 

It is now easy to see that the entanglement fidelity for a Gaussian channel has 
the form 

F{p,i) = J d 2 aG(a)\C p (a)\ 2 , (14) 

which is the same as the corresponding form IjlOJI of the pure-state fidelity. Notice, 
however, that the entanglement fidelity is not given by the mixed-state version of 
Eq. Q. A more thorough discussion of the entanglement fidelity for mixed input 
states is given in Ref. [T5] . 

Another possible fidelity measure for a mixed-state input regards the mixed 
state as coming from a particular ensemble of states, {p n , *& n }, by which we mean 
that in each experimental run, one of the states \^ n ) is selected randomly with 
probability p n . Thus, in each run, the Gaussian channel delivers with probability 
p n the initial state with fidelity J r (^ / n ,7). A statistical average leads to a mean 
ensemble fidelity given by 

F{{Pn^n},l) = Y,PnH*n,l) = J>„(* „|$ (| „><* „|) | „) . (15) 

n n 

The mean ensemble fidelity depends on the particular ensemble used to make up 
an input density operator p = ^2 n p n \^n){^n\- The convexity of the entanglement 
fidelity means that the entanglement fidelity for a mixed input state p is less than 
or equal to the mean ensemble fidelity for any ensemble corresponding to p. 
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5. Scaling law of channel fidelity 



If we assume that the input state is given by a pure state, the fidelity is given 
by the quantum overlap formula (JSJ), which leads to the equivalent phase-space 
formulas @ and ()10|) . Here we use these phase-space formulas to show that the 
fidelity of the Gaussian channel obeys a universal scaling law. 

Applying the Fourier transform relation © between the Wigner and Weyl 
functions to Eqs. © and ()1U|) . we obtain two new formulas for the fidelity. The 
result is four equivalent forms for the channel fidelity: 



The first two lines are rewrites of Eqs. © and (jlUj) . The third line comes from 
writing the fidelity (JHJ) as an overlap of the input and output Weyl characteristic 
functions; it is thus also obtained by Fourier transforming the Wigner functions 
in the integrand of the second line. Similarly, the last line is obtained by Fourier 
transforming the Weyl functions in the integrand of the first line. 

The first and third forms (and the second and fourth) show us that 



This scaling law reflects a duality between Wigner functions and Weyl charac- 
teristic functions. This duality is related to sub-Planck structures of the Wigner 
functions For a given input state, the Wigner function has two important 

scales: a small scale I and a large scale L. The small scale I characterizes the 
sub-Planck phase-space structures in the input state's Wigner function. The large 
scale L characterizes the scale over which the Wigner function is nonnegligible. For 
pure states, these two scales are related by an uncertainty relation, LI ~ 1. Since 
the Weyl characteristic function is the Fourier transform of the Wigner function, 
these two scales appear inversely in the Weyl function. 

The scaling law ((T7j) displays this duality. If we look at the fidelity in the form 
of the second relation in Eq. (|T6|) . we see that the fidelity between the input state 
and the output state approaches unity if the dispersion of the Gaussian channel 
satisfies a = -y/7/2 < I so that the integral approaches Tr(/9 2 ), which is 1 for a 
pure state; in this case, all the small-scale phase-space structure of the Wigner 




(16) 



^,7) = -^*, 4/ 7 ). 



(17) 
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function is well transmitted through the noisy channel. The dual form of the 
fidelity, expressed in the the last relation of Eq. ((16j). says that to get good fidelity, 
we need to have 1/a = -y/2/7 > L, so that the integral reduces to the square of 
the integral over the entire Wigner function. Putting these two results together 
gives the phase-space uncertainty relation, LI ~ 1. 

In terms of the dispersion, the scaling law has the form 

r(%*) = ^r(%i/a). (is) 

Additional discussion of the relation between fidelity and sub-Planck structure and 
of the phase-space uncertainty relation can be found in Ref. |15j . 



6. Maximum fidelity 

From these considerations, we can derive the maximum channel fidelity that 
can be achieved by any initial pure state. We show that the maximum fidelity is 
given by the coherent-state fidelity for all values of 7. To demonstrate this, return 
to the expression for the fidelity given by the last formula in Eq. (|16j) : 

^•(¥,7) = J d 2 ad 2 Pe-^ a -ft 2/2 W*(a)Wv(P). (19) 

The task can be restated as finding the pure state that maximizes this overlap. 

Notice that this fidelity can be thought of as the average value of e _7 ' a ~^ I 2 
with respect to a pure product copy state, \^ a)®\^b), of two modes, A and B; the 
joint Wigner function of the two modes is Wab(®, P) = W^{a)W^{(5). Introducing 
modes C and D, with annihilation operators c = (a + b)/\^2 and d = (a — b)j\pi 
and corresponding c-number variables x = ( a + P)/V2 and 5 = (a — /3)/v2, we 
can rewrite the fidelity as 

JF(* ;7 ) = J d 2 5e-^ 2 W D {5), (20) 

where Wd($) = J d 2 x Wab{x> What we see is that the fidelity is the expec- 
tation value of the mode-D operator whose symmetrically ordered associated 
function is e -7 ' 5 ' 2 . Letting 7 = (n + 1/2) -1 , we see that A^ is given by n + 1/2 
times the density operator for a thermal state of mode D whose mean number of 
photons is n = 7 _1 (1 — 7/2). Thus we can write the fidelity as 



^(^,7) =Tr{A lPD ), 



(21) 



[Author and title] 



8 



where 



1/2 + fi ( n ' '' 



1 + n V 1 + n . 

i ( i _ 7 / 2 y td 



1 + 7/2 Vl + 7/27 

1 /l -7/2\ ( at ^ fet )( a - 6 )/ 2 



1 + 7/2 VI + 7/2 



(22) 



Notice that Ay =0 = Id i confirming that the fidelity is 1 regardless of the input 
state for 7 = 0. Generally we can bound the fidelity by the largest eigenvalue of 

•7 r ( l I / ,7) < |max eigenvalue of Ay^ = - ^ — . (23) 

The reason this is the largest eigenvalue is that the factor in large parentheses in the 
expression for has magnitude < 1, which means that the largest eigenvalue, 
corresponding to the vacuum state for mode D, is (1 + 7/2)~ 1 . Since coherent 
states saturate the upper bound, we can write 

^(7) = ^- (24) 

The bound on the expectation value of A^ is useful in other applications than 
the noisy Gaussian channels considered here. For that purpose, note that the 
bound holds for all joint states pab of modes A and B — i.e., it holds for the 
expectation value Tt^A^pab) — not just for the pure product copy states that are 
relevant to the channel fidelity. A joint state achieves the expectation value bound 
if and only if mode D is in vacuum; i.e., the state pab is the state of the two 
output modes, a = (c + d)/y/2 and b = (c — d)/y/2, of a 50:50 beamsplitter that 
has vacuum incident on its mode-D input. 

If we specialize to pure product input states, \*$>a) <8> \&b), we can say much 
more about when the bound is achieved. Since mode D is in vacuum, we have 

= d\* A ) ® |*b> = ^(a - 6)|* A > ® |$ B > , (25) 

which implies that 

a\V A ) ® |$ B ) = \V A ) ® 6|$ B ). (26) 

This requires that a \ ^ a) = ($bH$b)|* a) = a I* a) and 6|$b) = (*a|o|vI/a)|$b) = 
a\^s}, i-e., that \^a) and |$b) are the same coherent state \a). Thus the maxi- 
mum channel fidelity is achieved if and only if the input state is a coherent state. 
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7. Examples of pure states in a Gaussian channel 

7.1. Number states 

Let us take as an example the case of an input state that is a number state 
|n) of a harmonic oscillator. In this case the fidelity of transmission through a 
Gaussian channel is 

^(|n),7) = — f d 2 ae- 2 ^ 2 ^\(n\D(a)\n}\ 2 . (27) 
Using the property 

(n\D(a)\n)=e-^ 2 / 2 L n (\a\ 2 ), (28) 

where L n denotes the nth-order Laguerre polynomial, one can calculate an exact 
expression for the fidelity generating function for all number states: 



^( 7 ,A) = X;A^(|n), 7 ) 



V(l + 7/2) 2 - 2A(1 + 7 2 /4) + W ~ ill) 2 

= 1 . (29) 

V[(l-A) + (l + A)7/2] 2 -A 7 2 

The resulting fidelity of a Gaussian channel with a number state at the input is 

(1-7/2)" p f l + J 2 /4 \ ,„ n * 

^ (|n) ' 7) = (i + 7/2)^ p " v r^ 2 74j ' (30) 

where P n {x) is a Legendre polynomial. For 7 = 1, 2, this becomes 

2 (2nV 
H\n)A) = ^+1^(5/3) , H\n), 2) = 22n \i[ nl) 2 ■ ( 31 ) 

Using a series expansion of the generating function ()29|) or working directly with 
the expression (|30j) . one can easily calculate the fidelities for the lowest number 
states: 

4 + ■y 2 16 + I67 2 + V 4 

^ 1 ^) = 2 (2^7)3' ^l 2 )^)=2 + . (32) 

Techniques similar to those used in this section can be used to calculate the 
channel fidelity for an input state that is an arbitrary superposition of number 
states, but in the absence of some general technique like the generating func- 
tion ((2*n|) , the calculation becomes increasingly tedious as higher number states are 
included in the superposition. As an example, for input state |^) = (|0) + |l))/v2, 
the channel fidelity is 

1 + 37/4 + 774 
^' 7)= (1 + 7/2)3 • (33) 
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7.2. Squeezed state 

A squeezed state of a one-dimensional harmonic oscillator is given by the fol- 
lowing formula |19| : 

| /i ) = (l-| M | 2 ) 1 /4 e -*«» t2 /2|o). (34) 

Here 

\fi\ = tanhr = \J ^ ^ _ , (35) 

where r is the usual squeeze parameter and n is the mean number of oscillator 
quanta. For \i = 0, n = 0, and the squeezed state reduces to the ground state of 
the one-dimensional harmonic oscillator. The fidelity of a Gaussian channel with 
a squeezed state at the input is 

.T(|/x),7)= ; 1 (36) 

KW sfl + (2n + 1) 7 + 7 2 /4 

Figure 1 depicts the fidelity of various states discussed in this section, as func- 
tions of the channel noise 7. 

8. Example of mixed-state fidelities for a Gaussian channel 

Because of the many experimentally uncontrollable properties of the input 
states to a quantum channel, in reality we have to deal with sources described by 
a mixed state or, in some cases, by a particular statistical ensemble of incoming 
pure states. In Sec. 14.2.1 we introduced the entanglement fidelity and the mean 
ensemble fidelity as fidelity measures to characterize these situations. We note 
that since the mean ensemble fidelity (|15[) is an average of pure-state fidelities, 
it satisfies the scaling law l[T7jl and the upper bound (|23|) . Moreover, since the 
entanglement fidelity is bounded above by the mean ensemble fidelity, we can 
write generally that 

Hp, 7) < ?({pn, * „}, 7) < (37) 

for any ensemble that corresponds to the input state p. 

In this section we illustrate the two mixed-state fidelities by considering an 
ensemble of number states, {pn,\ n )}> that are selected with the Bose-Einstein 
probabilities 

n n 

Pn = (l + ft)"+i • (38) 
The corresponding density operator is the thermal density operator 

p = y2p n \n)(n\ = — - ( — - ) , (39) 
^ l + n \l + nj 

n=0 
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Fig. 1: Plots of the channel fidelity as functions of 7 for different input states. 
From the top: the vacuum state |0), the squeezed state \fi) with n = 1, and the 
first two number states, |1) and |2). 



where n is the mean number of thermal quanta. 

In this case we can apply the formula (|29j) to find the mean ensemble fidelity. 
A simple calculation gives 

f = 1 = = 1 , (40) 

+ (2n + l) 7 /2] 2 - n(l + n) 7 2 y/l + (2n + 1) 7 + 7 2 /4 

the same as the fidelity (|36|) for an input squeezed state with the same mean 
number of quanta. It is easy to verify that this ensemble fidelity satisfies the scaling 
law (|17|). The entanglement fidelity follows from inserting the Weyl function for a 
thermal state, C p (a) = e~' a ' (n+i/2)^ j n ^ (|T1)) . which gives 

It is trivial to see in this case that J-(p, 7) < !F, with equality holding if and only 
if n = or 7 = 0. In Fig. 2 we have plotted, as functions of the channel noise 
7, the mean ensemble fidelity and the entanglement fidelity for an ensemble with 
Bose-Einstein statistics. 
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Fig. 2: Plots of the entanglement fidelity (lower curve) and the mean ensemble 
fidelity (middle curve) as functions of 7 for n = 1. For reference, the upper curve 
gives the fidelity of a coherent state. 



9. Physical examples of Gaussian channels 

9.1. Cloning 

As an initial example of a Gaussian channel, we now investigate a very simple 
experimental setup used for approximate cloning of the states of a field mode. The 
no-cloning theorem |20[l21[ l22] shows that a universal and faithful cloning machine, 
which would clone an arbitrary input quantum state perfectly, is incompatible with 
quantum mechanics. It is possible, however, to find imperfect cloning machines 
that copy quantum states with some loss of quantum fidelity. The simplest device 
designed to clone quantum states is a 50:50 beam splitter that is preceded by an 
amplifier with amplitude gain of \[2 to compensate for the reduction in signal 
amplitude at the beam splitter. This setup is depicted in Fig. 3. 

The input (pure) state of mode A is cloned into two outgoing clones in modes 
C and D. The setup of this cloning device requires only linear optics and linear 
amplification. As a result, the annihilation operators for the outgoing modes of 
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Fig. 3: Cloning experimental setup with an amplifier of amplitude gain y/2. BS 
denotes a 50:50 beam splitter. The annihilation operators of the input mode, the 
vacuum mode, and the modes that carry the clones are indicated. 



the clones are 



-^(v^a + «L P + v) = a + 7y=(4m P + v ) , 

-^=(V2a + a{ mp - v) =o+-^(aj Bnp -u), (42) 



where ai mp is a vacuum-noise creation operator that describes the noise introduced 
by the amplifier, and v is the annihilation operator for the vacuum-mode V that 
is incident on the unused port of the beam splitter. The whole dynamics of the 
cloning process, including amplification, is described by a unitary transformation 
in the extended Hilbert space Ha <8> Wamp <8> Ti-v of the input mode A, the amplifier- 
noise mode, and the vacuum mode V. As a result of this transformation, the 
reduced dynamics of such a cloning device corresponds to a Gaussian channel 
with 7 = 1. 

For all initial pure states, the cloning fidelity is given by 

•fcioncW = ^(*, 7 = 1) = , 4) , (43) 

where the latter equality is a consequence of the scaling law (|17[). The general 
fidelity bound ()23|) for a Gaussian channel gives us a bound on the cloning fidelity, 
-^clonc < •^ 7 max(7 = 1) = §, with equality achieved only for cloning of coherent 
states. 
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9.2. Continuous- variable teleportation 

As a second example, we consider a continuous- variable quantum teleportation 
protocol, which is effectively equivalent to a Gaussian channel [22|. Quantum 
teleportation is a process that can transfer an arbitrary quantum state from a 
system held by one party, usually called Alice, to a system held by a second party, 
usually called Bob 123 • The process requires a pair of systems, shared by 

Alice and Bob, in an entangled state — the entangled resource — and an amount of 
classical information transmitted from Alice to Bob. 

The entangled resource is ideally a pure two-mode squeezed state ^2] of two 
modes, A and B, which have annihilation operators a and b. Such states can 
be described by a Gaussian Wigner function, W A b(cx, 0), that is specified by the 
following nonvanishing second moments: 

n+i = Jd 2 ad 2 f3\a\ 2 W A B(a,f3) = Jd 2 ad 2 p\(3\ 2 W A ,B(a,(3), 

m = J d 2 ad 2 f3af3W AB (a,f3) = J d 2 ad 2 f3a*f3*W AB (a,{3). (44) 

The moments must satisfy n > and \Jn(n + 1) > \m\ to ensure that the Wigner 
distribution corresponds to a valid quantum state. The state is pure if and only if 
yj n(n + 1) = |m|, in which case it is a two-mode squeezed state. When the state 
is pure, the limit \m\ —* oo gives the original entangled state of Einstein, Podolsky, 
and Rosen [2S], with m negative leading to a Wigner function proportional to 
^(qa + Qb)S(pa ~Pb), which are the ideal correlations for the teleporation protocol 
we are considering, and m positive leading to a Wigner function proportional to 
fi(QA — Qb)$(pa + Pb)- The correlated state of Eq. (fH|) is separable (unentangled) 
if and only if n > \m\ |29| . 

The state to be teleported is an input pure state p of a mode V in Alice's 
possession, which has annihilation operator v. The protocol consists of (i) Alice's 
measuring the two (commuting) homodyne quadratures contained in the Hermitian 
real and imaginary parts of the operator v + a^; (ii) Alice's communicating to Bob 
the (complex) result £ of this measurement; and (iii) Bob's displacing the complex 
amplitude of his mode B by £. 

The efficacy of the protocol is quantified by the fidelity between the output state 
of mode B and the input state \ip), averaged over the possible measurement results. 
The teleportation protocol involves three modes in an extended Hilbert space 
TCv 8> 7~C A ®7~Lb- Reduction of this protocol to the Hilbert space of the incoming 
mode reduces the teleportation protocol to a Gaussian channel with quantum 
fidelity given by Eq. ©, where the Gaussian noise distribution is calculated from 
the following relation: 

G(v) = [ d 2 a d 2 [3 <5(/3 + a* - v)W AB {a, (3) . (45) 
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Simple calculation involving the Wigner function of the entangled resource leads to 
a Gaussian distribution for G{y), with the noise parameter given by 7 = 2[l + 2(n+ 
m)]. A separable resource has n > \m\ and, accordingly, 7 > 2. For y 7 n(n + 1) > 
m > n, the correlated state is entangled, but with the wrong sort of correlations 
for the protocol we are considering, so 7 > 2. For ^Jnin+l) > —m > n, the 
correlated state is entangled and < 7 < 2. Perfect teleportation is achieved if 
7 = 0, which corresponds to m = —\Jn{n + 1) — ► —00 (m + n = —1/2), i.e., to 
a pure entangled state with perfect EPR correlations. For pure input states, we 
conclude that the maximum teleportation fidelity for a given entangled resource 7 
is 

1 1 

1 + 7/2 2(1 +n + m) 

and is achieved if and only if the input state is a coherent state. 

We have used the fidelity bound (|2~4")) to get two other interesting results 
The first is that the maximum fidelity for teleporting a coherent state using the 
standard protocol, but with any separable state for modes A and B [not necessar- 
ily a state of the Gaussian form specified by the moments (|44j)]. is 1/2. The second 
result has to do with local hidden-variable models for continuous- variable telepor- 
tation. The teleportation of any Gaussian input state can be described within a 
local hidden-variable model, no matter what fidelity is achieved in the teleporta- 
tion; the hidden-variable model is based on the classical phase-space variables of 
the Wigner distribution. For non-Gaussian pure input states, we have shown that 
the value 7 = 1 plays a special role: each non-Gaussian pure input state |^) has 
its own threshold fidelity, ^"(^,1) < ^ r m ax(7 = 1) = 2/3, below which its tele- 
portation can be accommodated within an extended phase-space hidden-variable 
model and above which it cannot. 



10. Conclusion 

In this paper we follow Hall [H] in defining a Gaussian channel as one in which an 
input mode is subjected to random Gaussian displacements in phase space. Such 
channels arise naturally whenever field modes undergo linear optical transforma- 
tions, linear amplification, and measurements of quadrature components, provided 
that the overall channel preserves the mean complex amplitude of the input mode. 
We introduce the quantum fidelity for both pure and mixed inputs and derive the 
maximum fidelity that can be achieved over all input states. This bound — and 
related ones that might come from using similar theoretical techniques — should 
prove useful in analyzing the performance of Gaussian channels. 
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